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Abstract 

We study polynomials that are orthogonal with respect to the modified 
Laguerre weight z~''^'^'^e~^^{z—l)'^^ in the limit where n, N ^ oo with N/n — > 1 
and u is a fixed number in R \ Nq. With the effect of the factor {z — 1)^'', the 
local parametrix near the critical point z = 1 can be constructed in terms 
of ^'-functions associated with the Painleve IV equation. We show that the 
asymptotics of the recurrence coefficients of orthogonal polynomials can be 
described in terms of specified solution of the Painleve IV equation in the double 
scaling limit. Our method is based on the Deift/Zhou steepest decent analysis 
of the Riemann-Hilbert problem associated with orthogonal polynomials. 
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1 Introduction and statement of results 



1.1 Modified Laguerre weights 

Let a, b he real constants, and let iV be a positive number. Define the (complex- 
valued) weight function to be 

«;(2) = 2"e-^^(^-l)2^ (1.1) 

where z° and {z — 1)^^ are defined with cuts along [0, oo) and [l,oo), respectively. 
That is, = l^l^e^"^'^^^ with < arg^ < 27r, and {z - 1)^^ = \z - l\'^b^^ib^rg{z-i) ^^^^i 
< arg(z — 1) < 2tt. Let S be a contour in C \ [0, oo) that is symmetric with respect 
to the real axis, tends to infinity in the horizontal direction and remains bounded in 
the vertical direction. This contour divides the complex plane into two domains Q±] 
see Figure [H We consider monic polynomials 7r„ of degree n, 7r„(z) = 2;" + ■ ■ ■ , that 
are orthogonal with respect to w{z) on S, i.e., 

j nn{z)z^w{z)dz = 0, j = 0, l,...,n- 1. (1.2) 

This is an example of non-Hermitian orthogonality, for which there is no general 
existence and uniqueness result associated with such w{z). It will be part of our results 
that in the asymptotic regime that we will consider, the polynomial vr^ uniquely exists 
for n large enough. 

Imz A 




Figure 1: The contour S 
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For 6 = and a > — 1 (with a ^ Z) we may deform the contour of integration to 
[0, oo). Then the orthogonahty condition (11. 2p reduces to 

/"OO 

/ 7r„(a;)x^+"e-^^'cix = 0, j = 0, 1, . . . , n - 1, (1.3) 
Jo 

which shows that in this case 7r„ is related to the classical Laguerre polynomial with 
parameter a, 

7rn{z) = {-ir^A"\Nz). (1.4) 

For properties of the classical Laguerre polynomial lI^\z), see e.g. Abramowitz and 
Stegun P or Szego [2]. For 6 = and a < — 1, the orthogonality condition ( 11. 3p 
is not valid, but the relation (11.41) still holds true. Asymptotic properties of the 
Laguerre polynomials with large negative parameters were studied by Kuijlaars and 
McLaughlin [3l H] in the regime where N = n, n ^ oo, a — oo in such a way that 

hm = A 

exists. Then it was found that the zeros of 7r„ accumulate either on an open contour 
(if y4 > 1) or on a union of a closed contour and a real interval (if < A < 1); 
see Figure [2l For the special value A = 1, the zeros typically (but not always!) 
accumulate on the Szego curve S defined by 

S := {z eC: \ze^-^\ = 1 and \z\ < 1}. (1.5) 

In the special case that a = —n + v with ^ No = N U {0} fixed, it was also 
found E] that the local asymptotics of the polynomials 7r„ of (11.41) near the point 
z = 1 are given in terms of parabolic cylinder functions D^{z). The main tool in 
the analysis of |3l HI O [6] is the matrix valued Riemann-Hilbert (RH) problem for 
Laguerre polynomials and the Deift/Zhou method of steepest descent, which was 
introduced in [7] and first applied to orthogonal polynomials in [SI [9]; see also [10]. 
The parabolic cylinder functions appear via a construction of a local parametrix in a 
neighborhood of the point z = 1. 

The main effect of the extra factor {z — 1)^^ in the weight (II. ip is in the local 
behavior of the polynomials near the point z = 1. From the point of view of Riemann- 
Hilbert analysis it means that in the critical case A = 1 the construction of the local 
parametrix with parabolic cylinder functions will no longer work. It is the aim of 
this paper to show that in this situation the role of the parabolic cylinder functions 
is replaced by the \Ef functions associated with a special solution of the Painleve 
IV equation. To the best of our knowledge, this is the first time that Painleve IV 
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(1.3 
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Figure 2: Zeros of "^V^) ^oi n = 40, and A = 0.81 (left), A = 1.001 (middle), 
and A = 2 (right). 



functions are used in the construction of a local parametrix in a RH steepest descent 
analysis. 

We were guided in our approach by a number of recent advances in random matrix 
theory and the theory of orthogonal polynomials. In critical situations it was found 
that Painleve transcendents appear in the description of the local eigenvalue statistics 
of random matrices as well as in the local asymptotics of orthogonal polynomials 
near singular points. For example, the local eigenvalue statistics at the opening of 
a gap in the spectrum of unitary invariant random matrix ensembles is described 
by the Hastings-McLeod solution of the Painleve II equation pTl [T2| [T3] . The same 
function is well-known to play a role in the Tracy- Widom distributions for the largest 
eigenvalue of random matrices [H] and the problem of the length of the longest 
increasing subsequence of a random permutation [15]. Singular behavior at edge 
points may be described by special solutions of the Painleve I equation [16], [IT] , the 
second member of the Painleve I hierarchy [TH [19] , and the Painleve XXXIV equation 
[20] . The Riemann-Hilbert approach was used in all of these situations. 

In the statement of the results in this paper we will focus on the behavior of the 
recurrence coefficients in the three-term recurrence relation 

TTn+liz) = {Z- bn)TTn{z) - a„7r„_i(z) (1.6) 

in the asymptotic regime where n — > oo with 

a = -n + u, N = n + V2Ln^^^ (1.7) 

with b, V and L are real constants and independent of n. Then of course the recurrence 
coefficients depend on the constants 5, v and L, but we suppress it in our notation. In 
this regime it turns out that we see the appearance of the Painleve IV transcendent. 

As noted before, for 6 = the polynomials reduce to rescaled Laguerre polynomials 
(11. 4p . From the recurrence relation for Laguerre polynomials 

- zL^:\z) = (n + 1)4^1 W - (2n + a + \)L^:\z) + (n + a)L^:\{z) (1.8) 



4 



and (11.41) it follows that under the scaling (11. 7p the recurrence coefficients behave as 

a„ = - + 0(n"^/2), 6„ = 1 - + 0(n"^), as n ^ oo. (1.9) 

n n^l^ 

For general 6 G M and z/ ^ No, we prove that (except for a discrete set of values 
L) a„ and 6„ exist for n large enough, and satisfy 

, /,(!) 

an = — + 0(n-3/2), 6„ = 1 - + O(n-i), 

for certain explicit constants a*^^^ and fo*-^-*, which depend on 6, u and L, that are 
explicitly calculated in terms of functions associated with the Painleve IV equation. 
For 6 = they reduce to a'-^-' = v and fo*-^-' = \/2L. 

Before we can state our result we ffist discuss the Painleve IV equation. 

1.2 Painleve IV equation 

The Painleve IV (PIV) equation is defined as 

hr + ttM^ + W + 2(5^ + 1 - 2800)^ , (1.10) 



ds"^ 2u \ds J 2 u 

where and Goo are constants. It is known that all solutions u{s) of PIV are mero- 
morphic functions in the complex s-plane and all poles of u{s) are simple with the 
residue ±1; see [21], [221 123]. For more properties of the PIV transcendent, see a 
very good review article by Clarkson [23] and references therein. The ^P-functions 
associated with Painleve IV are solutions of the following system of linear differential 
equations (Lax pair) for \E'(/,s) given ffist by Jimbo and Miwa 



with 



9^ 5^ 

^-A^, -^B^, (1.11) 



A = (^l + s + j{e - K)^ a, + y (1 - ^) a+ 

+ - (i^-e- 600 + ^(^^-26) ) (7_, (1.12) 

y \ lu J 

B =las + ya+ + -{K -Q-Q^)a_, (1.13) 

y 
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where y = y{s) and K = K{s) are defined by 



and 



0-3 



y ds 
K 



1 
-1 



-u 



2s, 

du 
ds 



+ + 2su + 



1 
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;i.i4) 
;i.i5) 

;i.i6) 



The Painleve IV equation (11.101) is the compatibility condition for the system (11. lip . 
It should be noted that the Lax pair is not unique. For example, one can find another 
Lax pair in Kitaev [26] and Milne, Clarkson and Bassom [27j. 

To explain the special solutions of PIV that will appear in our results, we recall the 
RH problem for ^(/, s) associated with PIV. In the RH problem ^(A, s) is viewed as 



1 S2 

1 



1 

S3 1 

< 



1 

Sl 1 



1 S4 
1 



327ri6cx)0'3 



Figure 3: The contour and the jump matrices for the RH problem for ^' 



a function of A, with s appearing as a parameter in the asymptotic condition fll.23p . 
The constants Ooo and appear in the behaviors at infinity fll.23p and at the origin 
fll.24p . respectively. The jump matrices in the RH problem include four parameters 
(Stokes multipliers) Sj, j = 1, . . . , 4, satisfying 



1 + S2S3)e 



27ri0o 



[S1S4 + (1 + S3S4)(1 + SiS2)]e 



2cos27re. 



;i.i7) 



The RH problem for \E' is (see |28i Section 5.1] and f29]) 
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(a) s) is analytic for Z G C \ (M U i R); see Figure O 

(b) on the contours M U (with orientation as in Figure [3]), 



where 

(c) as / — s> cxd: 

where Z~®°° is defined with a branch cut along zM_; 

(d) as / 0: 

^(/,s)C"^r®"^ = 0(1) for Re / > and Iml < 0. (1.24) 

where is defined with a branch cut along i M_ and the connection matrix C 
is any invertible matrix satisfying 

5i525'3S'4e2"*®°°"3 = C-^e-^^'^'^'C. (1.25) 

Then, from Theorem 5.2 in [28], we obtain a solution of PIV (ll.lOp with constants 
9 and 9oo, by putting 

u{s) := -2s - ^log ((*-i)(s))^^. (1.26) 

The solution depends on the constants s/s in ( ]1.22p . which as explained should satisfy 
the equation (11.171) . The equation (11.171) is invariant under the transformation 

{si,S2,S3, S4} {dsi,d~^S2,ds3,d~^S4}, dj^O. (1-27) 
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M/_(/,s)52, 

$_(/,s)^4e2"*®°°"3^ 



for / G M+; 
for / G iR+; 
for / G M_; 
for / G iM_, 



;i.i8) 
;i.i9) 

;i.20) 

;i.2i) 
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1 

si 1 

1 

S3 1 



^2 
^4 



1 S2 

1 

1 S4 
1 



1.22) 



(/ + + + O (1)) e(4-0- /--3, (1.23) 



The \l/ function then transforms as \1/ i— c? '^a/^vl/c^'^a/s -vyj^ich leaves the solution (11.261) 
invariant under the transformation (11.271) . The connection matrix C transforms as 
C ^ d-^^^^Cd'"'-'/^. The PIV solution (fT:26|) is a meromorphic function of s and the 
RH problem for \E' has a solution if and only if s is not a pole of u. 

The special solution of (ll.lOp that will appear in our results corresponds to the 
special Stokes multipliers 

= _(g(eoo+e)7rj _ ^~{0oc+e)niy-0a^7Ti^ ^ _g(2eoo+e)7ri_ (1-28) 

where 

e = -b and = u + b. (1.29) 

It is easy to check that ( 11.28^ indeed satisfies the relation (11.17^ . Furthermore, by 
(ll.25p . the corresponding connection matrix C is 

and Ti and r2 are arbitrary constants satisfying r2 = 0,if6 = m/2,mGZ. Throughout 
this paper, we will choose ri = 1 and r2 = 0, which means 

C=f.g(2ei-ew (1-30) 



or after a transformation (I1.27P 



1.3 Statement of results 

The asymptotic formulas for the recurrence coefficients a„ and 6„ in (11.60 are given 
as follows: 

Theorem 1.1. Suppose 6 G R, G M \ No and put 

Q = -b, eoo = z^ + 6. (1.32) 

Let u{s) be the special solution of PIV corresponding to the Stokes multipliers (11.280 
and let K{s) be defined in U.15\) . Assume 

a = -n + u, N = n + V2Ln^/^ (1.33) 
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where L is not a pole of u{s). Then, for every large enough n, the monic polynomial 
TTn of degree n, satisfying 



7Tn{z)z^-''e-^'{z - If^dz = 0, J = 0, 1, . . . , n - 1, 



uniquely exists. The polynomials satisfy the recurrence relation (11.61) with recurrence 
coefficients a„ and bn satisfying 



„ = i(z/-ir(L))+0(n-3/2) 



as n ^ oo 



1.34) 



and, if K{L) ^ v. 



In 



-L 



+ 0{n 



as n ^ oo. 



1.35) 



In the Theorem ll.il f ll.35p still makes sense when u{L) = 0. In fact, from the PIV 
equation in f ll.lOp . it can be verified that 

u'{L) = ±46 = T^b if u{L) = 0. 



From ffTlSj) . then we have K{L) = or K{L) = 29 = -2b. Then by L'Hospital's 
rule, we get from (11.351) 



bn 

br,. 



1 + ^ 

:1 + ^ 



K'{L) _ 

' K'jL) 
2{iy+2b) 



L 



-L 



+ 0{n'^), if m(L) = and m'(L) = -46, (1.36) 
+ 0(n-i), if u(L) = and u'(L) = 46. (1.37) 



Note that the factor {u + 26) in (I1.37P can not be zero because of the assumption 
K{L) ^ z/ in ffL35D . 

When L is not a pole of u{s), it is still possible that K{L) = v. In this case, 
the RH problem for ^> in (I1.18p - (ll.24p is solvable. Moreover, for large enough n, the 
monic polynomial 7r„ still uniquely exists. However something strange occurs for the 
recurrence coefficient 6„. Note that (11.350 no longer holds when K{L) = v. In fact, 
as n oo, 6„ does not tend to 1 but to another constant, which means there is a 
sudden change from K{L) ^ v to K{L) = v. At this time, we can not explain why 
this phenomenon happens. 

In [3U] , Murata found the following Schlesinger transformation 



u [s] 



2K{s){K{s) + 2b 



u(s) K(s 



V 



(1.38) 



where u*{s) is a solution to the PIV equation with parameters 6 = —6, Goo = i^+b+1. 
Then, with f ll.38p . one can rewrite fll.SSp as 



whenever L is not a pole of u*{s). 

The case z/ G Nq is special and we do not understand it at this moment. That this 
is a particular situation can already be seen in the Laguerre case treated in [4J, since 
it corresponds to Laguerre polynomials with a zero of large order at the origin, and 
the RH problem stated in Section [2] below does not have a solution in case z/ G No 
and n > u, see Proposition 2.2 of [3]. 

Theorem 11.11 could still be true in case G Nq, but our method of proof fails. 
Indeed, in the transformation fl3.33p we make essential use of the fact that u is not 
an integer. Therefore we exclude the case z/ G No from our considerations. 

Using the Riemann-Hilbert method, it is possible to obtain asymptotic expansion 
for polynomials 7r„(2;) in all regions of the full complex z-plane as n — >^ oo. Near 
the point z = 1, the expansion involves the \Ef-functions given in fll.lip . We will not 
discuss it in this paper. 

For the polynomials, we restrict ourselves to the asymptotic zero distribution. 

Theorem 1.2. Under the same assumptions as in Theorem \l.l\ and K{L) ^ v, we 
have that all the zeros of polynomials 7^n{z) tend to the Szego curve S given in (11.50 . 
More precisely, for any neighborhood U{S) of S, there exists a positive integer hq such 
that for any n > uq, all zeros of nn{z) lie in U{S). 

Outline of the rest of the paper 

In Section [21 we formulate the RH problem for orthogonal polynomials 7r„, which is 
the starting point of the asymptotic analysis. We have to distinguish the two cases 

Case I: z/ > and z/ ^ N, 

Case II: z/ < 0, 

since the steepest descent analysis is different for the two cases. The two cases are 
dealt with in Sections [3] and H] respectively. Together they contain the proof of The- 
orem [LTJ In Section [5l we prove Theorem ll.2[ Finally in Section [6l we show that, 
with certain parameters 6 and 6oo, the special solutions of PIV studied in this paper 
are given in terms of parabolic cylinder functions. 




n 



oo, 



(1.39) 
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Since the two cases Case I and Case II are independent of each other, we are going 
to use the same notation for the functions and variables. We trust that this will not 
lead to any confusion. 



2 The RH problem for orthogonal polynomials 

Consider the following RH problem for a 2 x 2 matrix valued function F : C \ S 



(a) Y{z) is analytic for z G C \ S; see Figure [T] Here S is an oriented curve 
whose +-side and — side are on the left and right while traversing the contour, 
respectively. 



(b) Y{z) possesses continuous boundary values Y±{z) such that 

'1 w{zy 



Y+{z) = Y4z) 







for z E T,. 



(2.1) 



(c) Y{z) has the following behavior as z oo 



Y{z)={I + 0(- 



z"" 




By a by now standard argument, originally due to Fokas, Its, and Kitaev |l7j, the 
solution of the above RH problem, if it exists, is uniquely given by 



Y{z) 



n„,{z) 



Pn-l[Z) 



2m 



2m 



W[S)Pn-l[S} 

s — z 



ds 



(2.2) 



where 7r„ is a monic polynomial of degree n satisfying (11. 2p and Pn-i is a polynomial 
of degree < n — 1. The existence and uniqueness of the monic polynomial 7r„ satisfying 
( II. 2p is equivalent to the solvability of the RH problem. 

The recurrence coefficients and 6„ can be expressed in terms of the solution of 
the RH problem. 



Proposition 2.1. Assume that the RH problem has a solution Y and write 





Y{z) 



i+^ + ^ + o{z-': 

z z^ 



z" 
z- 



(2.3) 
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as z oo. Define 



and 



'n 




(2.4) 



(2.5) 



where the subscript ij denotes the entry of the corresponding matrix. If 7^ 0, 
then the monic orthogonal polynomials vr„+i, 7r„ and Hn-i (of degrees n + 1, n, and 
n — 1) with respect to w{z) on E, exist and they satisfy the recurrence relation 



Proof. The proof is standard, see [TO]. For the above precise form, see also [HI 



As aheady said before, we are going to apply the Deift/Zhou steepest descent 
method to the above RH problem. In Sections [3] we treat the case z/ > and in 
Section m the case u < 0. Before we embark on the steepest descent analysis let us 
say a few words about the method. The method consists of a number of explicit 
transformations, which in this paper take the form Y \^ U ^-^ T S ^ R, which 
lead to a "simple" RH problem for R. Since the works of Deift et al. [H [9l [10] it is 
now clear what the main steps should be in the analysis of the RH problem associated 
with orthogonal polynomials, namely 

• normalization of the RH problem with the so-called (^-function; 

• opening of the lens around the oscillatory region; 

• construction of global and local parametrices. 

In any extension of the method (such as the one in this paper) these issues appear 
and it is unfortunate that all details have to be checked each time, since we do not 
have, and maybe cannot expect, a master theorem giving general conditions under 
which the method will work. 

In the present paper we are dealing with orthogonality on a contour in the complex 
plane. To handle this feature we largely follow the steepest descent analysis done for 
Laguerre polynomials with large negative parameters in [3] (for z/ > 0) and in [3] (for 
z/ < 0). This will lead to an analysis on varying n-dependent curves depending on 
An = I — -. For n — 00 the curves tend to the Szego curve. The precise analysis 




(2.6) 



Proposition 2.4]. 



□ 
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is rather delicate near the critical point z = 1, since we have to arrive at a local RH 
problem that can be modelled by the RH problem associated with Painleve IV. A 
simpler and maybe more natural approach would be to avoid the varying curves and 
work with the Szego curve from the beginning. We tried to do this but we were not 
able to handle all difficulties this way. 



3 Case I: z/ > and z/ N 
3.1 Introduction 

We assume u > and u ^ N. Throughout we write 

t := — = 1 + V2 Ln~^^\ (3.1) 
n 

Let 

a I' , , 

An:= = 1 - -. 3.2 

n n 

Since we will concentrate on the asymptotics as n oo, then without loss of general- 
ity throughout this paper we assume n is large enough. Thus, An G (0, 1). Following 
the Riemann-Hilbert steepest descent analysis in [1], we define 



/3i,„ := 2 - A„ - 2v/l - A„ = 1 + - - 2./-, (3.3) 

n V n 



(32,n := 2 - A, + 2^1 - A„ = 1 + - + 2J- (3.4) 

n \j n 

and 

Rn{z) := sj{z-l3i,n){z-(32,n), Z E C\ /32,n], (3.5) 

where Rn{z) z as z ^ oo. Let 



(t^n{z) = \ r ^^ds, zeC\ ((-oo, 0] U [/?!,„, oo)), (3.6) 



where the path of integration from jSi^n to z lies entirely in the region C \ {{—oo, 0] U 
[Pi^n, oo)), except for the initial point The curves where Re0„(z) is constant are 
trajectories of the quadratic differential (see Strebel [31]) 

- ^ds^ = Jl^M^^^ds^ (3.7) 

s s 

which has two simple zeros at and /32,ra and a double pole at 0. In P] it is shown 
that there exists a unique curve ro,n as follows: 
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Definition 3.1. The contour ro,n is a simple closed curve encircling once, so that 

Re0n(z)=O for zeTo,n; (3.8) 

see Figure m 

The curve ro,n depends on n and it tends to the Szego curve S (\1.5h as n ^ oo. 

Lemma 3.2. As n oo, the curve T^ n tends to the Szego curve S. 

Proof. From (13.31) and (13. 4p . we know that and P2,n both tend to 1 as n ^ cxd. 
Then, by (13.51) . Rn{z) ^ z — 1. As a consequence, 

4>n{z) ^ - I ds = -{z — 1 — log z) as n — > cxo, (3.9) 

2 Ji s 2 

uniformly for z bounded away from and oo. Therefore, from the definition of F^n 
in Definition 13. H we know that Fo,„ tends to a simple closed curve Fo,oo encircling 
once, and 

Re(z- 1 -logz) =0 for z G Fo,oo- (3.10) 
From (II. Sp . we see that Fq^oo is the Szego curve S. □ 

In a probability measure was defined for each A < 1 supported on a set 
Fq U [Pi, P2], where Fq is a closed contour and [Pi, P2] an interval on the positive real 
line. In the present situation we have that A = A^ = 1 — ^ depends on n, but for 
each finite n we still start from the corresponding probability measure /i = /i„ which 
now depends on n. The measure is given by 

d^ln{y) = ^"y^^ '^^^^^ ^ 27rz ^"'^'^^^ X[/3i,n,fe,n] iy)dy- (3.ii) 

Using (13. 6p and (13.80 . one can indeed check that (13. lip defines a positive measure on 

^0,n U [/3l,n; P2,n]- 

3.2 First transformation Y U 

Expecting the zero distribution as shown in Figure [2] for < A < 1, we start by 
modifying the contour in the RH problem. Let 

E^:=Fo,„U[/5i,„,oo). (3.12) 

We assume (without loss of generality) that the contour S in the RH problem for Y 
was chosen so that Fo,n is contained in see Figures [1] and HI 



14 



Figure 4: Contours S and TF = Tq^^ U [Pi,n, oo) 




Figure 5: The contour S*^ = T^ n U [/5i,n, oo) for the RH problem for U 



Introduce U as 



{Y{z) for zen+u no, 

, fl wiz)\ 
Y{z) [ y] for 2 G \ ^]o; 

see Figure HI Then we obtain the following RH problem for U : 
(a) U{z) is analytic for z G C \ TF , see Figure 



(3.13) 
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(b) U (z) possesses continuous boundary values on satisfying 

/i ^-n+u -Nz( _ 1 \2b\ 

U4z) = U4z)r e ^^z i) ^ for^ero,„, (3.14) 

f/+(x) = f/_(x)Q for x G 1) U (1, oo); (3.15) 

(c) U{z) has the following behavior as z ^ oo 

(d) U{z) has the following behavior as z — 1: 

f/(^) Q 1)"'^ = 0(1) for z G C±, (3.16) 

where 

c„^^ = -e"^ and c,,_ = -e^"^"*"^. (3.17) 

In fl3.14p and fl3.16p . the factor {z — 1)^^ is defined with a branch cut along [1, oo) (as 
before). 

Because of the different choices of branches in the definition of w{z) we find that 

M;+(a;) -u;_(a;) = |xr"+"e-^"|s-lp^(l-e2'^"'), xe(/3i,„,l), (3.18) 
w;+(x) -t/7_(x) = |x|-"+'^e-^"|x-lp'(l-e2(^+''')"0, x e {l,oo). (3.19) 

Therefore, (13.151) takes on different forms on the two intervals (/^i,™, 1) and (1, oo). 

Condition (d) is crucial in the above RH problem for U. With conditions (a)-(c) 
only, we do not have a unique solution for our RH problem for U. To obtain a unique 
solution, we have to add the extra condition (d) that controls the local behavior 
as 2; 1. This local behavior (13.161) with the precise constants (13.171) will also be 
important in the construction of a local parametrix around z = 1 that will ultimately 
be our goal. 

3.3 The g and functions 

The measure /i„ in (13.111) gives rise to the so-called (7-function, which depends on n, 
9niz)= J log(^-s)rf/i„(s), ^eC\S^, (3.20) 

ro,nU[/3i,„,/32,„] 
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where the logarithm log(z — s) is defined with a branch cut along , going to the 
right from s to oo. Recalling the definition of (pni^) in (13.61) . we similarly define 



4>n{z) 



Note that 



V 



z e C \ (-00, /52,^ 



0„(2;) = 0n(-2) ± —T^i for z G C^; 
see O Lemma 3.4.2]. We also introduce 

gtA^) ■■= 9n{z) + {t- l)gn{z) 

and 

(ptA^) := 0n(^) + {t- := + (t - 

where 
and 



l{z-Rn{z)), ZeVloo 
}{Z + Rn{z)), Zeno, 



Rn[z) 



Z G C \ [/?!,„, /32,n]. 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



Note that g^^z) = 0(1/2;) as 2; — > 00. It is easily seen that the functions Qniz) and 
(pniz) satisfy the following properties. 



Proposition 3.3. (a) We have 

1 



2^ l+0°(^), zenoHC^. 



(3.27) 



(b) For z G , we have 



z G To,™ nc±, 



fc) We also have 



-2(I)14Z)=2(I)1_{Z), ZeiPl,n,f32,n]. 



Z, Z E ro,n U [Pl^n,P2,n], 

Z-2(t)l{z), Z G [/?2,n,00). 



(3.28) 



(3.29) 
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Then from the above proposition and Proposition 3.2 in [1], we get the following 
properties for gt^n, (pt,n, and 0t,„. 



Proposition 3.4. (a) There exists a constant It^n such that 

9t,n{z) = - {An log Z + tZ + It^n) + < 

2 [±iA„7rz + 0t,„(z), zenoHC^. 

(3.30) 

Here \ogz is defined with a cut along [0, oo). The constant It^n is explicitly given 
by 

k,n = '^9t,nM - {An log X„ + X^), 

where x„ is the intersection of To n with the negative real axis. 

(b) For z G T,^ , we have 

{-(f)t^n{z) - (pt,n{^)T'^h zGro,„nc±, 
-20t,„,+ (z) - 2Anni = 2(t)t,n,-{z) - 2A„7ri, Z G (/3l,n, /52,n], 
-27ri, z G [/52,„,oo). 

(3.31) 

(c) We also have 

{An log Z + tZ + It^n, Z G T^n, 

An\ogZ + tZ + Aniri + lt,n, ZG [f3l,n, f32,n], 

An log z + tz + Aniri + It {z), z G [/52,„, oo). 

(3.32) 

3.4 Second transformation U ^ T 

To normalize the RH problem at infinity, introduce the second transformation U ^ T 

as 

^3 



T(z) = ((-l)"(e-*"'^-e*"'')j e-2"'*'""^f/(z) 

^ g-n3t,„(^)a3g|ni,.„a3 (^(_ 1)« (g"*^^ - e'""^)^ (3.33) 



Note that it is important here that u is not an integer, since this assumption guaran- 
tees that e"""^^ — e'^^" ^ and therefore T is well-defined. 
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With the RH problem for U and Propositions 13.31 and I3.4[ direct calculation gives 
us a RH problem for T as follows. As before (2; — 1)^'' is defined with a cut along 
[1,00). 

(a) T{z) is analytic for 2; G C \ S*^; 
(b) 



■l)"e- 



(3.34) 



^+W = T_(.)r eUJwj' (3.35) 

^+W=^-W| ' ^2nCrW j> -e(l,/?2,n); (3.36) 

T+(^) = T_(^) ( J 1^ ), ze(/32,n,oo); (3.37) 

(c) T{z) has the following behavior as ^ — >■ oo 

T{z) =1 + 



(d) T{z) has the following behavior as 2; — 1: 

T(z) (^J 1)''^ = 0(1) for zeC^ (3.38) 

with 

= (-l)'^(e-''"* - e^"0~^ e2"^*-"-±(^)-"'*'"c„,±, (3.39) 
where c„^± is given in (13.171) . 

To derive the jump matrices for T in fl3.35l) - fl3.37p . one needs to make use of 
fl3.15p . (13.181) . (13.191) and Proposition 13. 4[ The conjugation with the constant matrix 

( (— l)"(e~*'^'^ — e™'') ) in (I3.33P may look a bit awkward since it makes the jump 
matrix in (I3.34p more complicated. We have introduced it in order to simplify the 
jump matrices (I3.35p -( l3.37p . We also want to point out that the entry ^^^^^^^^^'^ in 
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(I3.36P and fl3.37p will come out in the Stokes multipliers s/s in the RH problem for 
PIV later on. 

By using the relations between the 5^- functions and 0-functions in f l3.30p and fl3.27p , 
the constant Ct,± in f l3.39p can be rewritten as 

3.5 Third transformation T S 

From fl3.8l) . (I3.22p and fl3.24p . we know that, when |t — 1| is small, the real part of 
(pt,n{z) + (l)t,n{z) is close to for z G Fo^n- This means the jump matrix for T on Tq ^ in 
fl3.34p has oscillatory diagonal entries. To remove this oscillating jump, we introduce 
the third transformation and open the lens as follows. In a fixed neighborhood of 
2; = 1 let the curve Fi be defined as Im^; = arg^;, then we continue Fi such that 
it becomes a closed contour with the Szego curve S in its interior. Note that Fi 
is independent of n and t. Then Fi and Fq,™ divide the complex plane into three 
domains fio? ^1 and Q^; see Figure [61 Define 




Figure 6: Szego curve and the contour for the RH problem for S 



(-1)' 



(2-1)^ 



S{z) - T{z) I , ^y„+i eY^'-|7' (_i)ne-"('^t."W+r.n(2)) J ' ^ e ^0; (3.41) 

S{z) = T{z)\ e-"'-e-^' ^n(ch,.r,(z)+it.„.(z)) ij' Z & ^i] (3.42) 

(2-1)^" ^ V 

Siz) = T{z), z e ni. (3.43) 

Direct calculation shows that S+{z) = S^{z) for z G Fo,n- Therefore, S{z) has an 
analytic continuation across Fo,n and we obtain a RH problem for S on the contour 
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S'^ = Fi U [0, oo) shown in Figure [71 We use Qq to denote the bounded domain 
enclosed by Fi. 




Figure 7: The contour E"^ for the RH problem for S 

(a) S is analytic for 2; G C \ S"^, see Figure [7J 
(b) 



S^{z) = S^{z 

S4z) = S4z 

Sjf-{z) = S-{z 
S^{z) = S^{z 



1 







- fill ^L/TVZ 



,n{4>t,n{z) + 4>t,n{z)) ^ 



]_ [z — l)2bg-2n0t,„(2)^ 

1 

^2n<f>t,n,+ (z) _ ]^~j26 

g ^2n<j)t,n,^{z) j ' 

g ^2n<f>t,n-{z) 

2b sm{u+2b)n ^2bTri ^-2nit.n(z) ' 
1 



1 2b s'm{u+2b)TT ^2bTTi 
sin un 



1 \z - ^■"■"e 

I I sin uw 







z e Fi; (3.44) 

(0,A,n); (3.45) 

ze{Pi,n,l); (3.46) 

zeil,P2,n); (3.47) 

ze{P2,n,ooy, (3.48) 



(d) 



S{z) has the following behavior as 2; — 00 

Siz) = I + 0[^ 

S{z) has the following behavior as 2; ^ 1 

1 CtA^-l?' 



Siz) 



'1 ctA^-i?' 



S{z) 







0(1), z e l]^nC^, (3.49) 

0(1), zeninc^, (3.50) 
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where is given in (13.401) . 



The jump matrices in (I3.44p - (l3.48p look quite complicated. However, most of 
them tend to the identity matrix as n — > cxd at an exponential rate. To clarify this 
fact, we need the following proposition. 




Figure 8: S is the Szego curve. The dark region is the region where Re(z — 1 — log^;) > 
0, the white region is the region where Re{z — 1 — logz) < 0. 



Proposition 3.5. Let 6 be a fixed small constant and Bs := {z ^ C \ \z — 1\ < 6} . 
Then there exist r] > and e > such that for all t E M. with |t — 1| < rj and all n 
large enough, we have 

Re (f)t,n{^) < for z eViXBs 

and 

Re0t,„(z) > e{\z\ + 1) for z G (0, oo) \ Bs. 
The same inequalities hold for Re(/)t^n(^)- 

Proof. First recall from (13. 9p . that (f)n{z) — > |(-2 — 1 — log 2;) uniformly for z bounded 
away from and 00; see Figure [S] for the property of Re(2; — 1 — log 2). For z near 
and 00 we have 




-logz + 0(1) as z 0, 
z + 0(1) as z 00, 



both uniformly for n large enough. Since both the contour Fi and the constant S 
are independent of n, there exists an e > such that Re0„(z) < —e on Fi \ Bs and 
Re (f)n{z) > + 1) on (0, 00) \ Bs uniformly for n large enough. From the fact that 

0°(z) = — !lLJ. ag ^ _^ unifomly for all z, 
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and the definition of (pt,n{z) in fl3.24p . it follows by continuity that the same inequal- 
ities hold for (pt,n{z) if t is sufficiently close to 1. 

Since Re0„(2;) = Re0„(2;) for all z G C, the inequalities also hold if we replace 
0„(z) by (t)n{z)- □ 



3.6 Construction of the local parametrix 
3.6.1 RH problem for P 

From Proposition 13.51 we know that the jump matrices for S are exponentially close 
to the identity matrix if n is large, except for the ones in the neighborhood of z = 1. 
Now, let us focus on the RH problem of S restricted to a neighborhood Bs of z = 1. 
We seek a 2 x 2 matrix valued function P that satisfies the following RH problem. 

(a) P{z) is analytic for z & Bs \ S"^, and continuous on Bs \ S"^; 

(b) P{z) satisfies the same jump conditions on T,^ HBs as S does; see f l3.44p - fl3.48p : 

(c) on dBs, as oo 



b , 



P(z) = [I + [^] ] n^"^ uniformly for z edBs\ S^; (3.51) 



(d) P{z) satisfies the same local behavior near z = 1 a.s S does; see f l3.49l) and 
fl3.50p with the same constants Ct^± as given in fl3.40p . 

The factor n^'^^ in the matching condition fl3.51l) is special and unusual in the 
local parametrix construction. It comes out of our construction with the RH problem 
for \l/ from Section [L2] and we are not able to remove or simplify it. However, we can 
deal with the factor in the final transformation in Section 13.71 below. 

3.6.2 Reduction to constant jumps 

We look for P in the following form 

p( _ f^iQ(^)(^ie^"'^*^"^^^"^^^+^^^*^'''(^ - 1)"^''^^ for zeBsnn^, 

where ai ~ (^^^ ' defined with branch cuts along [1, oo) 

and (— oo, 1], respectively. Thus we have < arg(2; — 1) < 2tt, < arg(l — z) < 27r 
and 

arg(2; — 1) = arg(l — z) ^ 7ri for z G C^. (3.53) 
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The factors e-"'^*'"(^)'^3(z - l)''''^ and e-"^*-"(^)'"3(l - zf"^ in (l332|) are introduced 
to cancel the corresponding factors in the jump matrices in fl3.44l) - fl3.48l) and make 
them independent of z. The resuh is that Q will satisfy a RH problem with constant 
jumps. 

Reorient the contour Yp fl so that all parts are oriented away from the point 
z = \\ see Figure [91 Then it can be verified that, in order for P{z) to satisfy the 
desired RH problem, Q{z) should satisfy a RH problem as follows: 









1 \ 







Figure 9: The contour S"^ fl i?5 for the RH problem for Q 



(a) Q{z) is analytic for z G -B^ \ S*, and continuous on Bf,\Yr\ 



(b) 



Q+(^) = g_(z) 
g+(^) = Q_(^) 
Q+(^) = g_(z) 



1 0' 

sin(z/+2fe)7r ^(;^-|-fc)7ri ^ ' ' 



1 
1 0\ 



z G (1,1 + 5); (3.54) 
^ G Fl n ^5 n C+; (3.55) 
0G (1-5,1); (3.56) 



Q+{z) = Q^iz) r j e2(-+^)™3^ ^ G Fl n n C- 



(3.57) 



(c) for z G as n — > oo, 

Qiz) = (^/ + o(^-^^^ (v^(^-l))-(^+'')'^^exp 



n 



z — 1 — log 2;) 



+ - 1) + 2 + log'^ - 1) + 2 '^'^ ~ ^)^* 



0-3 (3.58) 
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where \ogz and (z — 1) ^'^^''^ are defined with cuts along (— oo,0] and along 
Fi n n C^, respectively; 

(d) Q{z) has the following behavior as z ^ 1: 

/ n (p-^T^i _ pV-ivi\pb-Ki\ 

Q{z)[ ^ ^ ^ Xiz-V)""^--^ = 0{\), ^G^]^^C+; (3.59) 



Q(^) I e-- ^ ^-2... |(^ - 1)''^^ = 0(1), z^^l^ C-; (3.60) 



pi^TTl p — l^TTZ 



Q(^) I ? 1 - 1)'"' = 0(1), ^ G n C+; (3.61) 



1 



g(^) I ,(2.+V. ° 1 (^ - 1)'"^ = 0(1), zen'^n (3.62) 



pfTTi p — i/TTl 



where now (z — 1)^ is defined with a cut along Fi fl i?5 fl C . 

In the above RH problem for Q, it is notable that, although P2,n G (1, 1 + ^), the 
jump matrix in (13.541) is the same for z < (32,n and z > f32,n- This is in contrast to 
the RH problem for P, where we have different expressions for z < /52,n and z > /32,n- 
A similar thing happens for the jump matrix in fl3.56p . 

To obtain the matching condition in (13.581) . we need to make use of the following 
asymptotic formulas for (p^{z) and (f)n{z) as n ^ cxd. For z G dBg, we have from 

= li^ - 1) - + 0(^] (3.63) 



2' ' 2n{z-l) 
and from (13. 6p . calculated with the assistance of Maple 



(Pn{z) = i(2: - 1 - log2;) - ^logn + - 21og(2: - 1) + 2ni + logz + logz/ - ij 



Z/2 



where log z and log(2; — 1) are defined with cuts along (— oo, 0] and [1, oo), respectively; 
see also [5] and [6]. For z edBsHQ^, we have from (Km . (133B and (13321) 

-ba3p{ni4>„i^) + it-l)<t>0„(z))-Uu+b)ni)a3 



Q(^)= (/ + 0(-^)) (v/^(^ - l))-^-e 
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Substituting (13.631) and (13.641) into the above formula, we get 



Q{z) =(^I + (v^(^ - exp 



ft 

-{z-1- log 2;) 



+ ^(2vr^ + log ^ + hgu - 1) + Oi-) + ^^^^ ({z - 1) + Oi-)) + b)m 

2 n 2 \ n J 2 

where log 2; and {z — are defined with cuts along (—00, 0] and Fi fl i?^ fl C~, 

respectively. Recalling the definition of t in (13.11) . the above formula immediately 
gives us (13.581) for z G dBs fl Qq. Similarly, we can prove that (13.581) also holds for 

z G OBs nn^. 

To obtain the limiting behaviors (I3.59p - (l3.62p near z = 1, some careful calcula- 
tions are needed. We take (I3.62p as an example and (I3.59p - (l3.6ip can be obtained in 
a similar way. From (I3.52p . (I3.53P and condition (d) for the RH problem for P, we 
have, as z ^ l,z & fl 

Q(z)e-("^*^"--«+5(-+^M-3((^_ *)fe-3 1^ ^ = 0(1), (3.65) 

\ct-{z-l) 1/ 

where {z — lY is defined with a cut along [1, 00) and q _ is given in (I3.40p . Changing 
the branch cut of (2 - l)Ho Fi n n C" we get 

(2 - 1) (z - 1)6^"* foTzen^nC-. (3.66) 

From (KE^ and (KE^ . it follows that, as ^ ^ 1, z e fi^ n 

g(^)e-M..,-(i)+|{-+''WV3((;,_l)e-)^-3(^^^ ^^=0(1). 

Right multiplying e^n<t>t,n,-ii)+liu+b)Tn)a3 ^-bmaa both sides of the above equation 
yields, as 2; — > 1, z G fl 

Q(z)(z-lf''-n , ^.,}, • \=o{\). 

' I Q _ — 1) e ''^ 1/ 

Recalling (I3.22p . (I3.40p and moving [z — Xf'^'^ to the right of the lower triangular 
matrix, we obtain (I3.62p . 
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3.6.3 Comparison with the RH problem for \1/ 

Now we compare the RH problem for Q with the RH problem for ^ in Section 11.21 
From fl3.54p - fl3.57p we see that we need the Stokes multipliers 

smz/TT (3.67) 

Because of the factors {z - if''^ in fl339|) - fl3:62|) and {z - l)-(^+'')'^3 f[338D . we are 
led to choose 

= -6, e^ = ;y + 5. (3.68) 

Then the above Stokes multipliers are obtained from those in f ll.28p . by applying the 
transformation in f ll.27p with 

Again it is important that u is not an integer, since this ensures d ^ oo. Furthermore, 
from fll.3ip . we get 



1 0\ / 1 



p(2i/ + 36)iri 



which appears exactly in fl3.62p . Thus the jump matrices and the local behavior near 
the point 1 in the RH problem for Q correspond exactly to the jump matrices and 
the behavior near the origin in the RH problem for 

From fl3.59p - fl3.62p . we do not want a logarithmic singularity at point 1 for Q for 
all values of b. This is special when we want to construct our parametrix for Q by 
using the -function in Section 11.21 Note that / = is a Fuchsian singular point of 

the equation — - = A'i! in fll.lip and the leading term of A as Z is 
ol 



7/7/ 2K 

(e - K)as - + —{K - 20)a_ 

2 uy 



The eigenvalues of the above coefficient matrix are B and — G. When = m/2, m G 
Z, usually this is a resonant case and there is a logarithmic singularity at for \E'; 
see [28] Section 1.3]. In this paper, we are going to make use of -function which is 
free from logarithmic singularities. Then with = —6, 0oo = u + b and the Stokes 
multipliers given in fl3.67p . in Section [H] we will see that there exist special function 
solutions to PIV when b = m/2, m G Z. Also see [32] for similar cases for Painleve H 
solutions. 
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3.6.4 Construction of Q 

Using tlie RH problem for we now construct Q as follows. We use the mapping 
function 

f{z) = [z-l-\ogzf' (3.69) 

which is a conformal map from a neighborhood of 1 onto a neighborhood of 0. We 
have 

f{z) = ^{z - 1) - y^{z - ir + 0{{z - If) asz^l. (3.70) 

Since for 2; G Fi fl -B^ we have that Im z = arg z, it can be shown that / maps Fi fl Bg 
to the imaginary axis. 
Now we define 

Q{z) = E{z)^ (^n'^fiz), L^^'j (3.71) 

where 

E{z) = i j {ue-^ ^)f<x3g^(.-b)™3^ (3 72) 

L is given in fll.33p and ^(-j ■) is the solution of the RH problem fll.18p - fll.24p with 
0, 9oo given in (13.680 and Stokes multipliers Sj given in fl3.67p . Because of fl3.70p we 
have that E is analytic in a (small enough) neighborhood of z = 1. We also see that 

z — 1 

L — ^ — ^ L cis z — ^ 1, 

V2f{z) 

Since L is (by assumption) not a pole of the PIV solution m(s), there is a small enough 
5 > (depending on L) so that 

z s = L- 



V2f{z) 

maps Bs to a compact subset of the s-plane that does not contain any poles of u{s). 
Then \1^(A,L ^j^^.^ ) exists for all z G Bs, and Q{z) is well-defined and analytic for 
zeBs\ s^. 

It is now easy to check that Q satisfies the jump conditions fl3.54p - fl3.57p . We 
also find that the behavior as z ^ 1 in fl3.62p is satisfied. This follows from the 
corresponding behavior as A in the RH problem for The behaviors as z — 1 
in fl3.59p - fl3.6ip follow from this and the jump conditions in the RH problem for \E'. 

So what remains is to check the matching condition fl3.58p in the RH problem for 
Q. For this we need to note that the asymptotic condition fll.23p holds uniformly for 
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s in compact sets away from the poles of u{s). Thus by the definition fl3.7ip we have 
for z G dBg as n — s> oo 



Q{z) = E{z) J + O ^ exp 



2/ W^3+ 



0^3 



(v+b)cFz 



where f{z) ^^^^'^ is defined with a cut along Fi n 5^ fl C . Since E[z) is independent 
of n, then 



Q{z) 



I + O 



1 



exp 



n 



n2 



■as 



1 



n^^fiz) 



{i/+b)a3 



Eiz). 



From the definition of E{z) in (13.721) . we get 



Q{z) = J + O 



{^{z - l))-('^+^)-3 



exp 



n 2. N n^\/2L(z- 1^ 
2^ ^""^ 2 



+ ^(log2; + logz/ - 1) + ^(z/ - h)^ 



(3.73) 



where log 2; and {z — 1)~('^+^) are defined with cuts along (—00, 0] and Fi fl -85 fl C~ 
respectively. Recalling the definition of f{z) in (13.691) . we obtain (13.581) from (13.731) . 
This completes the construction of Q. 



3.7 Final transformation 

Now let us consider the difference between the exact solution S and the parametrix 
P constructed. Define 

R{z) = { ^ ^ ' ^ ^ ' 3.74 

\m''-'S{z)P-\z), zeBs\J:^. 

From the RH problems for S and P, since P is analytic for z E Bs\T,^ and satisfies 
the same jump conditions on T,^ n Bs as S does, R{z) is analytic for z E Bs except 
for a possible pole at 2; = 1. Moreover, because 5* and P have the same local behavior 
near z = 1, R{z) has to be analytic at z = 1. Thus, we get a RH problem for R{z) 
on a contour as follows: 

(a) R{z) is analytic for z G C \ S^; see Figure [TOj 

(b) R+{z) = R4z)Jr{z) for z G S^; 

(c) R{z) = I + 0{l/z) as 2; 00; 
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Figure 10: Contour for the RH problem for R 



(d) R{z) = 0(1) as z^l. 
The jump matrices in the above RH problem are given by 

^^(^) = ( e-^^i-e-^^' -b„n(,i>t.„(z)+it.r,.(z)) i|' Z^TiXBs, 

= ( 1 J' ^ €(0,1-5), 

Q ' ' ^^"'^l 1, 2: e (1 + 5,00), 

J^(z) = P{z)n-i''\ zedBs. (3.75) 



From Proposition [33], the jump matrices are exponentially close to the identity matrix 
if n is large except for the one on dBs. For dBs, we have the following lemma: 

Lemma 3.6. The jump matrix Jr{z) on dBs is given by 

J^(z) = Piz^n-^"-' = 1 + -^P^-^\z) + -P'^-^^z) + 0(n-h, (3.76) 

m n 



uniformly for z E dBs, where P^ and P^ '^\z) are given by, for z E Bs 

p.-.,,) . [Lj^) ..^(., - (3.77) 

and 

^'^ I' Liz + I) ^ 1, , , fz-l\ , 

(2; -1)2 V2(2;- l)/(2;) V V2fiz)J 

(3.78) 
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Here E{z) is given by IjjS. 7S\ ), and \l/_2 0,1"^ given in U.23\) . 



Proof. Since Q{z) = E{z)'^ (n^f{z),L 



2-1 

'V2f(z) 



we have from ^M) and (l332i) 



Piz) = aiE{z)<iI in-^f{z),L 



z - 1 



(n(Mz)+{t-l)<P'i{z))-^{,^+b)7Ti)a3, 



z-1] 



-60-3 



(3.79) 

for z G BsHflQ. Using the formulas for 0°(-2) and 4>niz) in fl3.63l) - fl3.64l) . the definition 
of E{z) in (13.721) and the asymptotic formula for in (ll.23p . we get from (13.791) 



z-l 



X E{z) exp 



vL{z + l) 



0-3 + 



2n(^-l) n(z-l) 



nP{z) 
-(73 + 0(n- 



+ 0(n 



(3.80) 



which gives us (13.771) and (I3.78P for z E Bs\^VLq. Similarly, one can prove (13.771) and 
([3178]) also hold for ^ G 55 nfi^. □ 

For large n, the jump matrix is close to the identity matrix, both in and in 
L^-sense on S^. Then following similar analysis as in [9l [33], for large enough n, we 
know the RH problem for R is solvable. Moreover, from (I3.76|) . R{z) can be expanded 
as follows 

R{z) = 1 + -^R^~^\z) + -R^-^\z) + 0(r2-i), 



n 



n 



(3.81) 



as n ^ 00, uniformly for z G C \ S . To prove our Theorem II. 1[ we need to derive 
explicit asymptotic formulas for R^~^\z) and R^~ 



■2)1 



z] as z — 00. 



3.8 Asymptotic formulas for ^^(z) and ^^(z) 

To derive asymptotic formulas for R^~^\z) and R^~'^\z) as 2; — > 00, more information 
about \Ef_i(s) and \E'_2(s) in (ll.23|) is required. Using Proposition 1.1 in p. 51], 
one can derive the explicit formulas for them as follows: 



-H 



H 



(3.82) 



^-2(5) 



;^_2(s))ll (^-2(s))l2 
;^-2(s))21 (^'-2(S))22 



(3.83) 
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where 



(^-2 




1 

2 




2 2^ 




-i)(0 + eoo) 




»)l2 = 


y 

2 


(i— : 


). 








»)21 = 


K 

— {K-2Q)- 
uy 


y 


-e- 




(^-2 


»)22 = 


1 
2 


H^-sH- 


2 2^ 




-M)(e + eoo)" 



y and K are given in (11.141) and (11.151) . respectively, and 

H{s) =H := ^{K - 29) - + s) (ir - - Goo). 

Then, we obtain the following expansion for R^~^''{z) as 2; — > 00. 
Lemma 3.7. We have that 



(3.84) 
(3.85) 
(3.86) 
(3.87) 

(3.88) 



M 



M M 



H — - + 0{z ) as z ^ 00, 



z z^ 



where 



and 



M := ResP(-^)(2) 

z=l 



fH{L)-vL ^{K{L)-u)p- 
\-'^y{L)p^ -H{L) + uL 



(3.89) 

(3.90) 
(3.91) 



p = 2-('^+^)/2 (^e-i)<^/2 

Proof. From the RH problem for R and the expansions in (I3.76P and (I3.8ip , we know 
that R^~^\z) satisfies a RH problem as follows: 

(a) R^~^\z) is analytic for 2; G C \ dBs; 

(b) R^^^\z) = i?L"'^(z) + P(-i)(^) for z G dBs; 

(c) P(-^)(2) = 0(1/2) as 2 ^ 00. 

Here is clockwise oriented. From (13.771) and the property of f{z) in (13.701) . we 
know that P^~^\z) is analytic in except for a simple pole at 1. Thus, R^~^\z) is 
explicitly given by: 

^ -KesF''''{z), zeCXBs, 

(3.92) 



Z - 1 z=l 

^ Res P^-^Hz) 



p{-i)( 



zeC\Bs, 

z G Bs. 



Z -1 z=l 

This formula immediately gives us (13.891) . From (I3.77p . (13.820 and the definition of 
E{z) in fl3T2|) . we get (13:901) with p given by (IMTD . □ 
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Similarly, we get the expansion for '^Kz) as z oo. 
Lemma 3.8. We have that 

(z) = ^ + V = - + — —2 + 0{z-^) asz^oo, 3.93 

[z — 1)^ z — I z z^ 

where -B_i is a constant matrix and B^2 is given by 

^ ((V&-2(^))12 + f (^-2(i^))ll - ^^i^(^) + 

(3.94) 

p given by AS.Ql]) and (\E'_2(s))jj given in l^3.84\ )- (3^8l\ ). 

Proof. As in the proof of Lemma [3. 7[ we have a RH problem for R^~'^\z) from (13.761) 
and dMH): 

(a) i?(-2)(z) is analytic for ^ G C \ OBs; 

(b) /?t'^(z) = R^S^\z) + P^-^\z) + P(-i)(z)(i?(-i))_(z) for z G 95^; 

(c) R^-^\z) = 0(1 /z) as 2 oo. 

From (I3.70p and (I3.78p . we see that P'--^\z) is analytic in Bg except for a double 
pole at 1. If we rewrite P^~'^^{z) + P^~^\z){R^-^^)-{z) as 

p(-2)(^) + p(-i)(^)(i?(-i))_(^) = + ^ + W^(^), (3.95) 

[z — 1)^ z — 1 

where B_i and B_2 are two constant matrices, W{z) is an analytic function in Bs- 
Then, R^~'^\z) can be given as follows 

2; G C \ Pj, 



R(''){z) = { {z-iy z-V ' (3.96) 

-W{z), z G 



Thus, as 2; 00, we have 



R^-'\z) = + 0(^-^)^ + + ^ + ^(^"')) ^-1' (3.97) 

which gives us ( Km . From (l3T8l) . ([333]), (13:951) and the definition of E{z) in (l3T2|l . 
we have (Km . □ 
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3.9 Proof of Theorem 11.11 in Case I 

Then we are ready for the proof of Theorem 11.11 

Proof. Reversing the transformations Y U ^ T S ^ R, we can recover Y from 
R. Since the RH problem for R is solvable for large enough n, it follows that the RH 
problem Y has a solution for large enough n. This implies that the corresponding 
monic orthogonal polynomial 7r„(z) uniquely exists for large enough n. 

Next, we are going to calculate the asymptotic formulas for the recurrence coeffi- 
cients ttn and bn- From the definition of gn in (I3.20p . we have 

9n{z) = j \og{z - S)dfln{s) 

and so by (13.231) 

gt,ni^) = 9n{z) + {t- l)fi'°(2;) = log 2 - diZ~^ - rfg^"^ + 0{z~'^) as 2: ^ CX), 

for certain constants di and d2 (that depend on n and t). Then, from (13.331) . (13.431) 
and (I3.74p . we get 

R[z) = ({-ly^e-'^" - e''^'')n-')~'"'e-5"'*>"-^ + ^ + ^ + 0{z-'^)^ 

X (/ + + ^/ + + 0{z-')^ e^"''-- (^(_i)n(e-- - e-)n-') ^ 

= I+^ + ^ + 0{z-') (3.98) 
z z^ 

with 



= c-''''{Y_i + dina:^d'\ (3.99) 

R-2 = C"'^^(^y-2 + c/iny_iCr3 + ^1 + ^2^0-3 jc'^^ (3.100) 

where c = (^(-l)''(e-''''' - e^''^)^"^^ e^"'*.". Recalling ([23]) and ([23]), it immediately 
follows that 
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and 



12 



12 



12 



-R-1 ) + din 

22 



) + diTi 

22 



12 



22 



(3.102) 



12 



From (ESI]), (ESnD and (^M), we know that 







+ 




1 


i¥ 


n 





+ 



M 



0{n 



Combining (ICTl) and (I3l03i) . it follows that 



'n n 



and 



i?_2 = + - (e^i + + O(n-i). 

Thus, from (I3.10ip . we get 

an = -{M] [m] +0(n-5). 
n V / 12 V / 21 



(3.103) 

(3.104) 
(3.105) 

(3.106) 



From fl3.90p and the assumption that L is not a pole of u{s) and K{L) ^ u in Theorem 
rn we have (M)i2 = ^{K{L) - v)p~'^ ^ 0. Then, from fl3:T02|) . we have 



M 



12 



^((5_i)i2 + (5-2)12) +0(n-i) 



m) +-^(i?_0i2 + O(n-i) 



iM22 + ^(-"^ 



1 + 



12 

+ 



1 (5_i)i2 , 1 (B-2)l2 



M 



+ 0{n-^] 



1 + 



/n\ / 22 

1 (S-.)i. 1 /^A 

/n V / 22 



1 - 



1 



12 



n 



M 



+ 0{n-'] 



12 



n (M)i2 



(3.107) 
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With (K90\f and flSlMj) . the above formulas give us fOD and ffL35|) . 
This completes the proof of Theorem 11.11 for Case I. 



□ 



4 Case II: z/ < 
4.1 Introduction 

Since z/ < 0, with the definition of in (13.21) we know An > 1. Following the RH 
analysis in [3], define 

(3n:=2-An + 2%^ An - 1 (4.1) 

and 

Rn{z) := ^(z-/5„)(^-:^J, ^ G C \ ro,„, (4.2) 

where Rn{z) behaves like z as z —>■ oo. Recalling the value of An in (13. 2p . /5„ can be 
rewritten as 

Pn = l + - + 2m/ — . (4.3) 
n \ n 

As in [3], define 

</'n(^) := I r ^^ds, zeC\ (ro,n U Ti,, U [0, oo)), (4.4) 

where the path of integration from to z lies entirely in the region C \ (ro,n U Fi^^ U 
[0, oo)), except for the initial point Similarly define 

M^) -=1 i ^^ds, zeC\ (Fo,„ U F2,„ U [0, oo)), (4.5) 

where the path of integration from /5„ to z lies entirely in the region C \ (Fo,n U F2,n U 
[0, oo)), except for the initial point It is not difficult to verify (also see [3]) 

(l)n{z) = (j)n{z) =F TTZ for z G Vt^ . (4.6) 

To clarify contours Fo,n, Ti.n and F2,„ in (14.41) and (14. 5p . we consider trajectories 
of the quadratic differential 

- ^ds- = Jl^Ml^ds\ (4.7) 



which has two simple zeros at /?„ and /3„ and a double pole at 0; see Strebel [3T]. In 
[3J it is shown that there exist curves Fo,n, Fi^^ and F2,n as follows: 
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Definition 4.1. The contour Fq,™ is a curve from to f3n which crosses the negative 
real axis, so that 

Re0„,±(2) = O for ^ G ro,„. (4.8) 

Ti^n and r2,n are curves that form the analytic continuation of ro,n such that (j)n{z) 
and 4>n{z) are real and positive on them, respectively; see Figure fTTl 









r2,n 


ro,n X 


pr~~~~ — 




1 















Figure 11: Curves Fo,„, Fi „ and F2,n 

As in [3], in this case we choose the curve S in the RH problem for Y to be 
^ = ^o,n U Fi_„ U F2,n- And we do not need a step Y U as in the Case I. 
Following the similar analysis as in the proof of Lemma [3. 2[ we have: 

Lemma 4.2. As n ^ oo, the curve Fo,n tends to the Szego curve S U.5\) . 

As an analog of (13. lip in the Case I, we define a measure dfin to be 

dfiniy) = ^^^^dy, y G Fo,„. (4.9) 

Using (14. 4p and (14.80 . one can verify that d^n{y) is a probability measure on Fo,n; see 
also [3j. 

4.2 The g and functions 

As in the Case I, define the (7 -function to be 

gn{z):= [ log{z-s)dfin{s), Z G C \ (Fo,n U Fi,,) , (4.10) 

where the logarithm log(2 — s) is defined with a cut along Fo,„ U Fi Introduce 

gt,n{z) := gniz) + {t - l)gl{z) (4.11) 
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and 

:= 0„(z) + (t - l)0O(z), := + (t - 1)00(2;), (4.12) 

where 

gl{z):=\{z-R^{z)) and 0O(^) := ^ 6 C \ To,.. (4.13) 

Note that g^{z) = 0(1/2;) as z ^ oo. 

Using similar analysis as in Propositions 3.7 and 3.8 in [3], we get the following 
properties for gt^n of (j)t,n, which is an analog of Proposition 13.41 in the Case I. 

Proposition 4.3. (a) There exists a constant It^n such that 

gt,n{z) = ]^{An\ogz + tz + lt,n)-(t)t,n{z), z G C \ (10,^ U Fi,, U [0, CX)) ) . (4.14) 

Here log 2; is defined with a cut along [0,oo). And the constant U-n is explicitly 
given by 

kn = 2fi't,„(/3„) - (A„log/5„ + tpn)- 

(h) For z & T,, we have 

gt,n,+ {^) - 9t,n,-{z) = < , , , , ^ (4-15) 

(c) We also have 

{An log Z + tZ + k^n, Z G F^n, 

Anlog Z + tZ + lt,n-24)t^n{z), 2; G Fi,„, (4.16) 

AnlogZ + tZ + lt,n-'24>t,n{^)> ZeT2,n- 

4.3 First and second transformations Y T, T S 

To normalize the RH problem at infinity, introduce the first transformation Y T 
to be 

T(^z) = e~^"'*'"'^^y (;z) e""^*'"^^^'^^ g|n't,no-3 ^4 

From the RH problem for Y in Section [2] and Proposition 14. 3[ we obtain a RH problem 
for T as follows. 

(a) T{z) is analytic for 2; G C \ S, S = Fo,„ U Fi^„ U T2^n', 
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(b) 




]^ g-2n0t,„(2)/'^ _ j^\2f) 




for 2: G ri,„; 



for z e ro,„; 



(4.18) 



(4.19) 



1 



^-2ruf)t,„{z) ^ ^ _ -^\2b 



for z e r2,„; 



(4.20) 



1 



(c) as z ^ 00 



T{z) =1 + 



z 



1 



) 



zeC\^. 



From ( 14. Sp and ( ]4.12p . we know that, when |t — 1| is small, the boundary values 
of lie (j)t,n{z) on both sides of Fo^n is close to 0. This means the jump matrix for T 
on Fo,n in (14.181) has oscillatory diagonal entries. Similarly as in the previous case, to 
remove this oscillating jump, we introduce the second transformation and open the 
lens as follows. Let F4 „ be a smooth curve connecting /3„ and 1. Note that since 
Im Pn 7^ arg we can not choose F4 „ to be the curve Im 2; = arg z. However, as /3„ is 
close to 1 for n large, we can choose F4^ji such that it is close to the curve Im z = arg z 
for n large. Then choose T^^n be the continuation of F4^„ such that F3.„UF4^„ becomes 
a closed contour and contains Fo,n in its interior. Then the curves ^j,n,j = 0, 1, ■ ■ ■ ,4, 
partition the complex plane into domains floo and flk,k = 0,1,2,3; see Figure [T2l 
Define 




Figure 12: The contour for the RH problem for T 
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S{z) 



Tiz 



S{z) = T{z 
S{z) = T{z 
S{z) = T{z 



S{z) 



Tiz 



(2-1)2^ 

1 

1 

1 ^Z — l)26g-2n0t,„(2) 
1 



for z G f2oo 
for z G l^o; 

for 2; G fii; 



for z G 



for z E fl'^i 



(4.21) 
(4.22) 

(4.23) 
(4.24) 
(4.25) 



Direct calculation shows that S^{z) = S^{z) for z G ro,„uri^„ur2,n- Therefore, S{z) 
has an analytic continuation across ro,n U Fi^^ U r2,n and we obtain a RH problem 
for S on the contour S'^ = £3^^ U S4^„ U [0, 00). Here S'^ divides the complex plane 
into two domains Qq and fi^; see Figure [131 As before (z — 1)^* is defined with a cut 
along [1, cxd). 




r4 










1 1 



Figure 13: The contour T,^ for the RH problem for S 
(a) S* is analytic for z G C \ S'^, see Figure | 
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(b) 



^+W = ^-W((^_^)-2V,,4.) i) for.Gr3,„, (4.26) 

S+{z) = S.{z)\ e2"<^*.".+(^) e2"'^'.".-(^)) for (0,1), (4.27) 
,0 1 I 

^+W = 5-(^)((^_^)^26 ,-2n,,„(.)j forzGr4,„, (4.28) 



^+(^) = 5_(2)| e2"<^*-.+(^) e^"-/-'.-.-!^)) I forze(l,oo); (4.29) 
.0 1 



c as 2; — > cxo 



(d) as 2; ^ 1 



S{z) = 1 + 



S(')[l,_iy2, j = ^^^^ foTzeQ',nC^, (4.30) 

5(^)(^Q / ^ j = 0(1) ioTzeninc^, (4.31) 

S{z)r / M = 0(1) for^Gl^^nC-. (4.32) 

Like in the Case I, condition (d) is important for the above RH problem for S. It 
makes sure that we have a unique and desired solution for S. 



4.4 Definition of i/j functions 

For future analysis, it is convenient to introduce new branch cuts for the functions 
{z — 1)^^, Rn{z), 4't,n{z) and 0^ ^(2;). More precisely, we want the branch cuts taken 
on T^^n- Let Tj^^ and Fj _ denote Tj fl C"^ and Fj fl C~, respectively. 

First for {z — 1)^*, instead of taking the cut along [1, 00), we take new cut along 
^4:,n,- U F3 „ _ U {—00, Xq], where Xq is the intersection point of T^ n and (—00, 0]; see 
Figure UM Thus, 

■ v~ forzGfi^nC", 
{z-1)^ {] [ °° (4.33) 

elsewhere. 
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^^^^^^^^^^^^^^ 




r4,- 


r3.- 





Figure 14: The new cut of (2; — 1)^* 

Then, define 

Rr,{z) := y/(z-/3„)(z-^„), ^ G C \ 14,., (4.34) 
where Rn{z) behaves hke z a.s z ^ 00. Also introduce 



i 


> 


r4 






r4 




V 


1 




1 



Figure 15: The cuts of ilJt,n and ^/'t^ri 

V-t.nl^) := M^) + 1)^°(^) and := i>^{z) + (t - 1)V^°(;2) (4.35) 
where 

M^)-=l[ ^^ds, 2GC\(r3.n,-Ur4.„U(-oo,0]), (4.36) 

^n(^) / ^^ds, zec\ (r3.„,+ u r4.„ u (-00, 0]) (4.37) 
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and 



(4.38) 



For illustration of the branch cuts of ipt,n and ipt,n, see Figure [151 Here are some 
useful properties of these new auxiliary functions, which can be obtained directly 
from Proposition 4.6.1 in [5]. 



Proposition 4.4. (a) The connection between Rn and Rn: 



Rn{z) 



Rn{z), forzeQo, 
elsewhere. 



Rn{z), 



(4.39) 



(b) The connection between il)t,n and (pt,'. 



<Pt,n{z), 

-0i,n(2;), 

-'Pt,n{z) + v4„7rz, 

(t)t,n{z) - An-Kl, 
(t)t,n{z) + A„7rii, 

'Pt,n{z) + (2 - An)^, for 2; e r^oo n C" 



for zeiVl-iUVliU fioo) n C+, 
forz e nc+, 
for z e^lof] C~, 
for z G 

for z eViin C~, 



(4.40) 



(c) The connection between ipt,n and ipt,n-' 



i't.niz) 



■ipt,n{^) - Tl'^h for Z G Q^, 
Aniz) + z:Tri, for z G (7^. 



(4.41) 



Then, with the help of fl4.33p and Proposition 14.41 the RH problem for S can be 
rewritten as follows. Now (2; — 1)^^ is defined with a cut along r4^„ .UFs^^ _U(— cxd, xq]. 



(a) S is analytic for 2; G C \ S"^; see Figure [T3l 
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(b) 



S^{z) = S4z) ( _ J) , ^ e r3,.,^-, (4.42) 

S^{z) = I _ , ze rs,n,-, (4.43) 

■ 2nV't,n,+ (2;) n \ 

S4z) = 5-(^) (^(^ _ j , - e r4,,+, (4.44) 



^+(-) = ^-(-) ( _ ) ' - e r4,„,_, (4.45) 



S^{z) = S.{z) ^ ^ ^ ^ ^j, ^e(0,l), (4.46) 

5+(^) = 5_(^)( ^ ' M, ^e(i,oo), 

(4.47) 

(c) as z — > oo 

5(^) = / + O 

(d) as 2; — > 1 

^(^) _ l)-2. j ^ ^^^^ for ^ e Qo^ n C^ (4.48) 
Siz)i^^ / ^ j = 0(1) forzeQ^nC+,(4.49) 

J ' 1 j ^ ^(^^ forzeQ^nC-.(4.50) 

Following similar analysis as in the Case I, one can see that most of the jump 
matrices in the above RH problem tend to the identity matrix as n — > 00 at an 
exponential rate. More precisely, we have the following proposition. 

Proposition 4.5. Let 6 be a fixed small constant and := {z E C \ \z — 1\ < 6}. 
Then there exist rj > and e > such that for all t with |t — 1| < and for all 
n large enough, we have 

Re'il)t,n{z) < -e for z e {r^^n U r4,„) \ Bg 
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and 

Rei)t,n{z) > e{\z\ + 1) for z e (0, oo) \ Bg. 
The same inequalities hold for He tlJt,n{z) ■ 

Proof. This proposition is similar with Proposition 13.51 and the proof is similar, too. 

□ 

4.5 Construction of the local parametrix 

From Proposition 14.51 we know that the jump matrices for S are exponentially close 
to the identity matrix as n — > oo, except for the ones in the neighborhood of 1. Let 
us focus on the RH problem of S restricted to a neighborhood Bs of z = 1. We seek 
a 2 X 2 matrix valued function P that satisfies the following RH problem. 

(a) P{z) is analytic for z ^ Bs \ T,^ , and continuous on Bs \ S"^. 

(b) P{z) satisfies the same jump conditions on S'^fli^^ as S does; see ( I4.42p - fl4.47p . 

(c) on dBs, as n — > oo 

p(z) = ( I + O { ^)) n^"' for z e dBs\J:^. (4.51) 



(d) P{z) satisfies the same local behavior near 1 as 5* does; see f l4.48p - fl4.50p . 

As in the previous case, we can not avoid the factor raa""^ j,^ f |4.5ip . 

To transform the RH problem for P into a RH problem for Q with constant jump 
matrices, we seek P in the following form 

faig(z)(Tie("'^*-(^)-^('^+'')^^)'^3(z - l)-'"^» for zE BgOn^, 
P(z) = < - -I (4.52) 

[aig(^)(Tie("^''"(")+2('^-^)-^)-3(z - l)-''-^ for z e Bs H Q^, 

where {z — 1)^ is defined with a cut along Ti^n- U Ps,?!,- U (— oo,xo]. The factors 
^ni^tA^>3(^z - ly""'^ and e"'^*'"(^)'^3(2 - l)-^''^ in fH32|) are introduced to cancel the 
corresponding factors in the jump matrices in (14. 420 - 04. 470 and make them indepen- 
dent of z. 

Reorient the curves S*^ fl Bs and let them all extend from the point 1. Then it 
can be verified that, in order for P{z) to satisfy the desired RH problem, Q{z) should 
satisfy a RH problem as follows: 

(a) Q{z) is analytic for z E Bs\ S'^, and continuous on Bs \ S'^; 
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(b) 



Q+{z) = Q_{z) 
Q^{z) = Q_{z) 
Q+{z) = Q_{z) 
Q^{z) = Q_{z) 



1 

1 

1 

^gi^TTi g-i/7ri^g-(2!/+b)7ri 

g(3i^+fe)7ri 

1 



ze(l,l + <5); (4.53) 

2; G (r3,„,+ u r4,„,+) n Bs; (4.54) 

J, (1-5,1); (4.55) 

2(.+6)™3^ ^ ^ (r3_^^_ u r^_^^_) n b,- (4.56) 



(c) for 2; G dBj as n — i> cx), 

Q(^) = (/ + o(-^))(v^(^-i)) 

n3v/2L(2- 1) 



exp 



n 



-(2; -1 -log 2) 



z/ h 
-{\ogz + log(-z/) - 1) - -vri 



0-3 



(4.57) 



where logz and {z — 1) '^'^^^'^ are defined with cuts along (— oo,0] and r4^„_ U 
r3^„_ U (—00, xo], respectively; 



(d) Q{z) has the following behavior as 2; — > 1: 



Q{z) 










(u+b)iTi 







(-.- 


(1^+6)714 


Q{z) 






Qiz] 





1 

g(i/+6)7ri 



1 

{u+3b)Tri 



z — 




= 0(1), 


zen-^r]C+ 


z — 




= 0(1), 


z G ^0 n c~ 


z — 


J^^feo-3 


= 0(1), 


zeninc+ 


z — 


]^^feo-3 


= 0(1), 


zeninc- 



(4.58) 
(4.59) 
(4.60) 
(4.61) 



where {z — 1) is defined with a cut along _ U Fa „ _ U (—00, xq]. 



The above RH problem for Q is very similar to the RH problem in Section 13.6.21 
One can see that there is only a constant difference between corresponding formulas 
in these two RH problems. Also again there is no difference for the jump matrices 
( I4.54P and (I4.56P in the neighborhood of j3n and respectively. 
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To obtain the limiting behavior as z ^ 1 in fl4.58p - fl4.6ip . one have to make use 
of (14. 6 p and Proposition 14. 4[ To derive the matching condition in fl4.57p . we need the 
following asymptotic formulas for ip^i^) ^^^d ipn{z) as n ^ oo. For z G S-B^, we have 
from fOHl) 

and from f l4.36p . (calculated with Maple) 

^i,n(^) = ]^{z-l- \ogz) - ^log^ + ^( ~ ^^°S(^ -1) + ni + \ogz + log(-z/) - 1 
+ -,7^ + 0(^], (4.63) 



n'^{z — 1)^ 

where logz and log(z — 1) are defined with cuts along (— oo,0] and r4 „_ U r3^„ _ U 
(— oo,Xo], respectively; see [Sj. The rest of the calculations is the same as in the Case 
I and we do not go into details. 

Now we are in a situation very similar to the Case I in Section [3.6.31 First, from 
the factors {z - l)^'"^ in fH38|) - fl4:6T|) and {z - l)-('^+^)<^3 fQTD . again we have 

e = -b, e^ = u + b, (4.64) 

which are the same as in fl3.68p . Next, from fl4.53l) - fl4.56l) we need the Stokes multi- 
phers 

_ -{u+b)m _ {u+3b)m _ _„{u+b)iTi 

' *2 - e , 

— ^gi^i"* Q-uwi^^-{2u+b)ni _ ^(3u+b)TTi \ " / 

Although (13.671) and (14.651) are not the same, the above Stokes multipliers, like those 
in fl3.67p . are again equivalent to those in (11.281) under the transformation in fll.27p . 
but with a different constant d = —e~"'^^. Then following almost the same analysis 
as in the construction of Q in fl3.7ip . except for a constant difference, we construct 
our parametrix by 

Q{z) = E{z)^ (n-2f{z), L^^) (4.66) 



V2f{z) 

with the same f{z) as given in fl3.69p and 

E{z)=lj^j (-z/e-^ z)t-3e-2™^ (4.67) 

There is one thing remaining about the curve (Fa „ur4 „)n-B5 under the mapping of 
f{z). Observe that f{z) does not map (F3 „UF4 ,i)n-B5 to the imaginary axis. However, 
for n large enough, according to our choice of F3„ and F4„, /((F3_„ U F4„) fl Bs) 
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is close to the imaginary axis. In the RH problem for \E' in Section 11.21 by the 
principle of analytic continuation, we can modify i M to a new contour A„ such that 
/((Fs „ur4 „)nS5) C A„ for n large enough. Therefore, we can still use the parametrix 
constructed in fl4.66p . 

4.6 Proof of Theorem 11.11 in Case II 

Having Q as given in (14.66^ . we can continue the analysis as in Sections l3.7ti3.9i The 
rest of the proof is the same as in the Case I and we do not go into details. 



5 Proof of Theorem 1.2 



Then we are ready to prove Theorem II. 2[ 

Proof. Let us consider Case I first. From (12.21) . fl3.13p and fl3.33p . we have 

7r„(z) = (r(z))n = (f/(^))ii = (T(z))ne"^'-(^). (5.1) 
First for z G we get from (I3.43P 

7r„(^) = (5(^))ne"^*'"(^\ ^ G fi^. (5.2) 
By fl3T4D and ([Ml]), it is known that 

S{z) = n~i''m{z)ni'''' = rT^^"'' (/ + o(^)) n^"^ uniformly for ^ G C \ 5^. 



(5.3) 

Combining (15.20 and (15. 3p gives us 

7r„(z) = ( 1 + O ( ^ ) ) e"^'*'"(') uniformly for ^ G \ B^. (5.4) 



So, for n large enough, there is no zero of vr„(2;) for z G fi^ \ B^. 
Then let us consider z G fig • From (I3.22p and (13.410 we have 



(T(z))n = (-l)"e-2"'^'-(^)^^-^(5(z))n + (-l)"^^^3Y^(5(z))i2 for z G fioHC^, 

(5.5) 

and from flX^ and ^^m^ we get 



(T(^))n = (5(^))ii + e^"'^'-"(-)^-"(g(^))i2 for ;2 G fii n C±. (5.6) 
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Then, combining flS.ip . (15. 3p . (15. 5p and (15.60 gives us 



n''{z - 1) 



2b 



iR{z)) 



12 



for zeQonC^X Bs, (5.7) 



and 



^ngt,„{z) 



^^^^ '^j + nf'{z-iyf> 



„~U-Kl _ VTTl / 1 



for 2 G ^]l n \ ^5. (5.^ 



For {R{z))i2 in (EZD, we get from (ESI]), (ICTll and (13:921) 



(i?(z)) 



12 



V^(z - l)y(L) 



(^fs:(L) - z/) ^ 0(^-1) for 2 e \ Bs. (5.9) 



According to our assumption K{L) ^ u in Theorem 11.21 we know that 
{R{z))i2 7^ foT z E Qq\ Bs and n large enough. 



(5.10) 



Let U{S) be the neighborhood of the Szego curve S. Note that Fo.n tends to S as 
n — > oo; see Lemma [3.21 Using the similar analysis as in the proof of Proposition 13.51 
it can be shown that there exists e* > such that Re0j,„,(-2) > £* for z E Qo\ U{S). 
This means that e"^"^'^* "'^^) is exponentially small as n — >^ oo. Then, from (15.71) we 
have 



,(z) = e''^^'"(")(-l)"- 



z-1 



\2b 



V2(k{L) 



V 



^{z - l)y{L) 



p-^ + 0{n-^) 



(5.11) 



for z eVLq\ {Bs U U{S)).. So, for n is large enough, there is no zero of 7r„(-2) in this 
region. 

For z E U{S), similarly it can be shown that e"'^*'"^^-' is exponentially small 
as n — > oo. Thus, we get from (15.80 



n 



zEn,\{Bsuu{S)), 



which means there is no zero of vr„(z). Thus, we can see that all zeros of tt„{z) 
accumulate in Bs U U{S) as n — oo. 

For Case II, the analysis is similar. This completes the proof of Theorem 11.21 □ 
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6 The solution to PIV with special parameters 



It is a well-known fact that for certain parameters G and Goo, the PIV equation has 
special solutions given in terms of parabohc cylinder functions; see e.g. [22l [30l [311 
[35l [36l [37] . The special solution of PIV that is of interest in this paper (that is, the 
one characterized by the Stokes multipliers (11.281) ) turns out to be of this type, for 
certain values of b. We will discuss it briefly in this section. 

We recall the RH problem for from (I1.18l) - (ll.24l) with parameters 

G = -6 and Q^ = u + b, (6.1) 

and Stokes multipliers fll.28p given in terms of b and u by 

^ ^-{2u+b)ni _ ^-bni g^ ^ _g{2!/+fe)7ri K^-^) 

It turns out that for b G Z/2, the corresponding special solution of PIV can be 
expressed in terms of parabolic cylinder functions. If z/ G N, the parabolic cylinder 
function reduces to a Hermite function and the PIV solution is a rational function, 
see also [30l [HI [36l [3H] . 

We discuss the cases 6 = and 6=1/2 here. 



6.1 The case 6 = 

Denote by vl>(^='=/2) ^j^e solution of the RH problem for \E' in Section [L2] with param- 
eters G and Goo given in (16.11) . Stokes multipliers given in (16.21) and 6 = |. 
When 6 = 0, by ([6ll]) and 1^ we have 

G = 0, Goo = (6.3) 

and 

= e^^^* - 1, S2 = 1, S3 = - 1, S4 = -e^"^'. (6.4) 

Substituting ([63]) and ([63]) into ([rT8D - f[T:24D . we get the RH problem for ^'(''=0). For 
this special case 6 = 0, the RH problem is solved explicitly in terms of parabolic cylin- 
der functions -Di/(C), where Dy{Q is the solution of the second order linear differential 
equation 

/(c) + - ^ - \^") y(0 = (6.5) 

uniquely characterized by the asymptotic property 

^.(C) = CeK^ (l - + ^ (^)) C - oo, I argCI < ^; (6.6) 
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see [H Chapter 19]. Other solutions of the differential equation (16.51) are D^{—Q and 
D_^_i{±i(). These solutions are related as follows 



DJ-z) 



{e--^^'D^,^i{iz) + e^^'D_,_i{-iz)) 



D,{z) = e^^'D.i-z) 



e^^'D_,_,{-tz) 



(6.7) 
(6.8) 
(6.9) 
(6.10) 



see also [39l p. 117]. When i/ ^ Z, the solution to the RH problem for '^^^ is given 
as follows, which can be verified with (l6.6l) -( l6TT0l) . 



Re / > 0, 



£3 

2-'^e"^(l-e-2'^"0 ) V(/)(l-e-2^"0^, Re/<0, 



(6.11) 



where 



e-t-D_,(-zy2(/ + s)) -^D,^i{V2{l + s)) 
^i^e-t-D_,_i(-zV2(/ + s)) D,iV2il + s)) 



for argZ G (0, | 



e^2^W_,{-tV2{l + s)) ^D,_^{-V2{1 + s)) 
^e^-^D_,_i(-^v^(/ + s)) D,{-^{1 + s)) 



V{1) = 



e-'^-'D^S^il + s)) ^D,_i(-v^(/ + s)) 



^^e-t-D_._i(zV2(/ + s)) D,{-V2{1 + s)) 



e^^'D^,{iV2{l + s)) 
^e§-D_._i(?v^(/ + s)) D,{V2{1 + s)) 



2n 



for arg/ G (|, vr), 



for arg / G (vr, ■^), 



for arg / G (f ,27r); 



see also [5i i6j|. Then, by (16.111) and (I6.12p . we have 
(v^(^=°)(/))i, = -2te-^^ 



r(z/) (e2^^^ - 1) 
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(6.12) 

D^^i{V2{l + s)) for Re / > 0. (6.13) 



Recalling fl3.82p and (16.61) . we get from above formula 



y 



i-*oo r(z/) (e2^^* - 1)' 



= coe-'\ (6.14) 

where 



= r(.)(e^---i) - ^^-'^^ 



Furthermore, by (I1.14p we get 



m(^=o)(s) = -2s - ^log?/(''=°)(s) = 0, (6.16) 

which is the trivial solution of the PIV equation with parameter = 0. 

It is easily seen that, with (16.161) . the expressions (11.341) and (11.351) for the recur- 
rence coefficients of the generalized Laguerre polynomials indeed reduce to ( II. 9p . 

6.2 The case b = l 

Knowing the precise form of \Ef^^^°^(/), we are going to compute \E'*^^"^/^''(/). For & = |, 
by (16.11) and (16. 2p we have 

e = -^, eoo = iy + l (6.17) 

and 

s^ = -ie^"''' + i, S2 = i, S3 = -ie"'^"''' + i, 3^ = -ie^"^'. (6.18) 

Define 

$(/) := e-3™3^(^=l/2)(;)e3-^-3(^(6=0)^;))-l_ ^g_ig) 

From the RH problems for vl>(*=°) and vl>(*=^/^), with l^/"^ defined with a cut along 
iM__, we easily obtain the RH problem for /^/^$(/) as follows: 

(a) is analytic for / e C \ (M U iM); see Figure El 

(b) {l^'^^l))^ = for leRUtR; 

(c) as / — s> oo 

where \E'_i(s) is given in (13.820 : 
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(d) as / ^ 

l^/^^l) = 0(1). (6.21) 

From parts (b) and (d) of the RH problem it follows that is an entire function 

in the complex A plane. By ( 16.201) we then have for / G C, 

Recalling (^L^"^^^^)i2 = -\y^^=^/^'^ (see (ESS])), we obtain 

lV^^{l)= I' 2y , (6.23) 

where * denotes entries that are not important for what follows. By fl6.19p . we have 
^1/2^(6=1/2) =e3™«/^/2<l>(/)^(^=°)(/)e-^"*"3_ (g_24) 

Combining fIgTTD . (ICTj) and (lO^ . we get that the (1,2) entry oil^''^^!^^=^/^\l) for 
Z in the fourth quadrant is equal to 



2^e-^-^ ^^^^^y^_^^ J._i(V2(/ + s)) - y-^^2'-^e-^^'DA^2{l + s)). (6.25) 
From fll.24p . we know 

^1/2^(6=1/2) ^ Q Q 1^ as Z ^ and / in the fourth quadrant. (6.26) 
This, together with (16.251) . gives us 

2'^-*"rM(fe)«-'(^^' - ^^2-lel.^Z,„(v/2.) ^ 0. (6,27) 

SO that 

From fll.l4p . we then see that the relevant solution of PIV in case 6 = 1/2 is equal to 

d 
ds 



«(fe=i/2)(,) = _2s--l logy (^=1/2) (s) 



d , D,iV2s) 



(6.29) 



log 



ds D,.i{V2sy 
From this formula, it is readily seen that the poles of 

^i(fe=i/2)(5) are the zeros of 

D^{y/2s) and D^^i{V2s). For u is real, the parabolic cylinder function 0,^ has 
max([z/],0) zeros on the real axis, where [■] is the ceiling function; see |39i p. 126]. 
So, the PIV solution ti(^=^/2) has poles on the real axis if z/ > 0. 
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6.3 The case 6 = |, A; G Z 

Repeating the previous calculations, it is possible (at least in principle) to find the 
solution ii(''='=/2) (s) for PIV with A; e Z in terms of parabolic cylinder functions. For 
example, for 6 = 1 we have 

^ ^ w(d,_,{V2s),D,{V2s)) 
/'='\s) = Cl— (6.30) 



«2 



where W is the Wronskian with respect to s and 



^ -r(.)(e---i)- ^^-^^^ 



Thus, 



^ >v(d,(V2sO,D.+i(v^s)) 



= ;7- log ^ — ^- (6.32) 

>v(l',_i(V2s),L',(V2s)) 

Again we see that there are poles on the real axis ii u > 0. 
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